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Abstract. ^ 

We point out that the amplitudes for the decays B — > D^D~ and Z? s — > D + D have no factor- 
izable contributions. If one or two of the D-mesons in the final state are vectors (i.e D* 's) there 
are relatively small factorizable contributions through the annihilation mechanism. The dominant 
contributions to the decay amplitudes arise from chiral loop contributions and l/N c suppressed 
\ tree level. We predict that the branching ratios for the processes BP d — > D^D~, B® — > D^~*D~ and 

l/l ■ B° d -> D+D7* are all of order (2 - 3) x 10 -4 , while 5° -> £>+£>^, 5? -► £>J*£> rf and D^D^* 

are of order (4 - 7) x 1(T 3 . If both D -mesons in the final state are D* 's, we obtain branching ratios 
of order two times bigger. 
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^ INTRODUCTION 

a: 

Qh| Decay modes like B — > mi, Kit are intensively studied. From the theoretical side, for 



instance within QCD-factorization. The decays of a B- meson into two D-mesons are 
different because the energy release is only of order 1 GeV, and therefore (QCD-) 
factorization is not expected to hold. B — > DD decays are also studied experimentally 

£0]. Here we discuss non-factorizable contributions to the decay modes B®, — > D^Ds*\ 

where D y s ' is a pseudoscalar or a vector meson. At quark level such decays occur 
through the annihilation mechanism bq — > cc, where q = d,s respectively. However, 
in the factorized limit the annihilation mechanism will give a vanishing amplitude due 
to current conservation (similar to Z) u — > K°K° Jl), unless one or two of the D-mesons 
in the final state are vectors. The contributions due to the annihilation mechanism are 
proportional to a numerically non-favored Wilson coefficient. In contrast, the typical 
factorized decay modes which proceed through the spectator mechanism, 5° — > D + Dj 
say, are proportional to the numerically favored Wilson coefficient. 

In our approach yl |4|], the non-factorizable contributions are coming from chiral 
loops and from tree level amplitudes generated by soft gluon emision forming a gluon 
condensate. The gluon condensate contributions can be calculated within a recently 
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developed Heavy Light Chiral Quark Model (HL^QM) [5]. This model has been applied 
to processes involving 5-mesons in Both the chiral loop contributions and the 

gluon condensate contributions are proportional to the numerically favorable Wilson 
coefficient. 



FRAMEWORK 
Effective Lagrangian at quark level 

The relevant effective Lagrangian at quark level reads: 

J*f w = -^ Vcb V c *Xai(n) Q,(n) , (1) 

where q = d,s and a,-(ji) are Wilson coefficients that carry all information of the short 
distance physics above the renormalization scale /!. The matrix elements of Qi(fJ-) con- 
tain all non-perturbative, long distance physics below fi. Within Heavy Quark Effective 
Theory (HQEFT) the effective non-leptonic Lagrangian Jzfvv can be evolved down to the 
scale /i = A x ~ 1 GeV ji . 

The numerically relevant operators in our case are 

Q\=4{q L 'fb L ){c L y il c L ) , Q 2 = A{c L fb L ) (q L y^c L ) , (2) 

where L denotes a left-handed particle. At /i = Ay, which by construction is the 
matching scale within our approach one finds a\ ~ —0.35 — 0.07/ and a 2 ~ 

1.29 + 0.08/. Note that the Wilson coefficients a, are complex below \x = m c [8]. In the 
next subsection we will see how the currents in the operators in © are bosonized. 

In order to obtain all matrix elements of the Lagrangian dU) we need the Fierz 
transformed version of the operators in ©. To find these, we use the relation: 



+ (3) 



where i, j, I and n are color indices running from 1 to 3 and t a denotes the color matrices, 
a being the color octet index. One obtains 

2f = T^<22 + e 2 , Ql = ^rQx + Qu (4) 

ISc 

where the superscript F means "Fierzed", and 

Q x =A{q L ft a b L ) {clJ^cl) , Qi = 4{c L ft a b L ){q L Y ll t a c L ). (5) 



These operators generate contributions proportional to the gluon condensate. 



Heavy Light Chiral Perturbation Theory 



The Heavy Quark Effective Theory (HQEFT) Lagrangian is: 



^HQEFT = ±Qt ) iv-DQt ) + 



[m 



Q 



(6) 



where (x) is a (reduced) heavy quark field (b or c in our case) with velocity v, and 

Qv \x) is the field of a heavy anti-quark (c in our case). Furthermore, mg is the heavy 
quark mass, and is the covariant derivative containing the gluon field. 

After integrating out the heavy and light quarks, the effective Lagrangian up to 
0(mQ ) can be written as a kinetic term plus a term describing the chiral interaction 
between heavy and light mesons 



St* Tr 



(7) 



(±) 



where H a is the heavy meson field containing a spin zero and a spin one boson: 



(8) 



Here a,b are flavor indices and P± = (1 ±y - v)/2 are projecting operators. The axial 
vector field sdu in is defined as: 



(9) 



where ^ = exp[i(U/f)]. Moreover, / is the bare pion coupling and U is a 3 by 3 matrix 
which contains the Goldstone bosons 7l,K,r} in the standard way, and is the axial 
chiral coupling. 

Based on the symmetry of HQEFT, we obtain the bosonized currents. For a decay of 
the bq system (see Fig. El left) we have o]: 



Tr 



(10) 



(+) 

Vb 



where (up to QCD and 1/mg corrections! 5]) a# = fii\JmH for// = B,D. Further, Q 
is the heavy &-quark field, is its velocity, and fit is the corresponding heavy meson 
field. For the W -boson materializing to a D, the bosonized current qZT^ Qv e ^ i s a ^ so 
given by (TTOb . but with //i replaced by ' representing the D meson. Qv t ^ is the 
field of the heavy c quark, and Vc is its velocity (see Fig.[H right). 
The bosonized b — > c transition current in Fig.[T]is given by 



(+) 



<(C0)Tr 



H 



(11) 
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FIGURE 1. Factorized contribution for B° — > D + D s through the spectator mechanism, which 
does not exist for decay mode B° — > DfDj. There are similar diagrams with vector mesons. 




FIGURE 2. Factorized contribution for B° — > DfD s through the annihilation mechanism, 
which give zero contributions if both Df and D~ are pseudoscalars. 

where £(co) is the Isgur-Wise function for the B — »■ D - transition, and v c is the velocity 
of the heavy c-quark. Furthermore, a = vt, ■ v c = v/, ■ Vg = Mb/(2Md). For the weak 
current for DD production (corresponding to the factorizable annihilation mechanism 

in Fig. El), the current Lq[, £ ^ is given by (fTTT) with W + ' replaced by \ and 

C(co) is replaced by £(— X), where A = v 5 ■ v c = [M|/(2M^) - 1]. Note that is a 

complex function which is less known than £(co). 

The factorized contributions for the spectator and annihilation diagrams are shown in 
the Figs. [T] and I2J The first diagram does not give any (direct) contributions to the class 
of processes we consider, but is still important because it is the basis of our chiral loops. 

The chiral loop amplitudes visualized in Fig.|3]are of order (g^mK/^nf) 2 compared 
to typical factorizable amplitudes in processes where these exist. For instance , the 
ratio between the chiral loop amplitude for 5° — > DfDj and the factorized amplitude 
for 5° — > D + Dj is ~ —0.20 + 0.26/ (before the difference in KM structure is taken 
into account). For vectors (i.e. D* 's) in the final state there are similar diagrams with 
various combinations of pseudoscalars and vectors in the loop, but the diagrams in Fig.|3] 
constitute the two classes of diagrams 0]. 

In a complete analysis, counterterms to the chiral loops has to be included. These 
counterterms are not considered here (or in 0,|j]) and has to be be considered together 
with the constant (non-logarithmic) chiral loop terms which we also have dropped in this 
analysis. The inclusion of counterterms and constant chiral loops terms will be discussed 
elsewhere. 



FIGURE 3. Non-factorizable chiral loops for B° — > ZJ+D S . There are similar diagrams for vector 
mesons in the final state. The two diagrams illustrate the two classes of diagrams. 



One might also write down possible terms consistent with HQEFT and chiral symme- 
try, for instance the following three terms: 



Tr 



^o^ a LH [ h 



■Tr 



m +) y a LH [ d V 



Tr 
■Tr 



Tr 



H^fLH^Yx 



(12) 



Such terms do not appear in the factorized limit, and will correspond to (at least) 1 /N c 
suppressed terms. Within pure chiral perturbation theory their coefficients are unknown, 
but they might be calculated within the HL^QM, as described in the next subsection. 



The Heavy Light Chiral Quark Model 

The HL^QM Lagrangian is 

■^HL^QM = HQEFT + ^%QM + ^Int ■ (13) 

The first term is given in © and the second term is described by the Chiral Quark Model 
of the light sector JToll involving interactions between quarks and (Goldstone) mesons: 

= X [f-(i@n + Y5^n) -m]z- (14) 

Here m — (230 ± 20)MeV is the SU (3) invariant constituent light quark mass, and % is 
the flavor rotated quark field given by Xl = £ 1l an d Xr = where q T = (u,d,s) is 
the light quark field. The covariant derivative contains the soft gluon field forming 
the gluon condensates (besides some chiral interactions) ifsL I^L FToh . 

The interaction between heavy meson fields and quarks is described by ^ : 



- „(±) n (±) , n (±)rr(±) v 



(15) 



where the coupling constant Gh = y/2mp/f, and p is a hadronic parameter depending 
on m (numerically p is of order one). For further details, see ref. ||5J]. 

The gluon condensate amplitudes can be written, within the framework presented in 
the previous section, in a quasi-factorized way as a product of matrix elements of colored 




FIGURE 4. Non-factorizable contribution for B° — ► DfD~ through the annihilation mechanism with 
additional soft gluon emision. The wavy lines represent soft gluons ending in vacuum to make gluon 
condensates. 



currents in ©, as visualized in Fig.lU The left part of Fig.|4]corresponds to the bosonized 
colored current: 



'■z?*aj*n(+)\ Gh8s G a 



Vb ho 64k » v 



xTr 



SYLH^W -F{o» v ,y. Vb }) 



(16) 



where G^ v is the octet gluon tensor, and F = 1%f 2 / (m 2 N c ) is a dimensionless quantity 
of the order 1/3. The symbol { , } denotes the anti-commutator. For the creation of a DD 

pair in the right part of Fig. HJ the colored current ^Qvt^ t a L()i s is bosonized 

similarly to ( fToT ). but involves He an d H^ . 

Multiplying the two colored currents, and introducing gluon condensate contributions 
by the replacement: 

oc 1 

*?G5vC%, - 4n 2 (^G 2 )-(g^ agv p -gtfgya) , (17) 

we obtain a bosonized effective Lagrangian which is \/N c suppressed compared to 
the factorized contributions. This effective Lagrangian corresponds to a certain linear 
combination of a priori possible \/N c suppressed terms at tree level (in the chiral 
perturbation theory sense). Among these are the three terms in (fT2l) . 

I /mi, suppressed terms seem to be negligible. In order to include \/m c terms, one 
must for instance consider the cc production current to this order: 

A7 M (cc) = — Q^ ) (ir- < D±(v c )fL+fLiY-D A _(vc))Q[;\ (18) 

where D^(v) = (g MV — v tl v v )D v . Within the HL^QM one may estimate both factor- 
izable and non-factorizable \/m c corrections due to this current There are also other 
operators. Compared to other contributions studied here, the relative size of l/m c sup- 
pressed contributions will be of order m/m c , where m is some hadronic parameter within 
the HL^QM with dimension mass, such as linear commbinations of m and (qq) / f 2 . The 
total contributions from such terms are significant and will be studied elsewhere. 



RESULTS AND DISCUSSION 



In our calculation we used the following input parameters: ccb — CCd — 0.33 GeV 3 / 2 , 
G H = 7.5 GeV" 1 / 2 and (f G 2 ) = [(315 ±20) MeV] 4 iBH, = 0.6, /„ = 93 MeV. 
We find ||4j] the following branching ratios 

Br(B°^D+DJ) =2.5 x 10~ 4 , Br{B Q s — > D + D ) = 4.5 x 10~ 3 , (19) 

Br(B° ^D+*D^) = 3.3x 10~ 4 , flr(fl° -> ) = 6.8 x 10~ 3 , (20) 

Br(B° ^D+D^*) = 2.0x 10~ 4 , 5r(5° -> D + D *) = 4.3 x 10~ 3 , (21) 

Br(B° ^D* + DJ*) = 5.4x 10~ 4 , flr(fl° -> D* + D *) = 9.1 x 10~ 3 . (22) 

The contribution of the constant term and the corresponding counterterm can change 
the branchig ratio for 5-meson decaying into two pseudoscalars by about 10%, while 
in the case of decay into one pseudoscalar and one vector D-meson, this contribution 
is in the range of 20 — 40%. In the case of 5-meson decaying into two vector mesons, 
the constant term is estimated to be 2-8 times larger than the logaritmic contribution, 
depending on the choice of the scheme in which the products of two Levi-Civita terms 
are considered. The uncertainty in input parameters can result in an additional error 
for the branching ratios. We estimate that it can be of the order of 20%. Within our 
approach the 1/mg corrections, with Q = c,b have so far been omitted. At least the 
l/m c corrections will be numerically significant. 



ACKNOWLEDGMENTS 

J.O. Eeg is supported in part by the Norwegian research council and by the European 
Union RTN network, Contract No. HPRN-CT-2002-003 1 1 (EURIDICE). S. Fajfer is 
supported in part by the Ministry of Education, Science and Sport of the Republic of 
Slovenia. 



REFERENCES 

1. For a recent study, see K. Abe et al., e-Print Archive: hep-ex/0508040 

2. J. O. Eeg, S. Fajfer, J. Zupan,Pfrys. Rev. D 64, 034010 (2001). 

3. J.O. Eeg, S. Fajfer, A. Hiorth, Phys. Lett. 6 570(2003)46-52. 

4. J.O. Eeg, S. Fajfer, A. Prapotnik, Eur.Phys.J. C42 (2005) 29-36. 

5. A. Hiorth and J. O. Eeg, Phys. Rev. D 66 (2002) 074001, and references therein. 

6. A. Hiorth and J. O. Eeg, Eur.Phys.J. direct C30 (2003) 006. 

7. J.O. Eeg, A. Hiorth,and A.D. Polosa Phys.Rev. D65 (2002) 054030. 
J.O. Eeg, K. Kumericki, and I. Picek, e-Print Archive: hep-ph/0506152 

8. B. Grinstein, W. Kilian, T. Mannel, and M.B. Wise, Nucl. Phys. B363 (1991) 19. 
R. Fleischer, Nucl. Phys. B 412 (1994) 201. 

9. R. Casalbuoni, A. Deandrea, N. Di Bartelomeo, R. Gatto, F. Feruglio and G. Nardulli, 
Phys. Rep. 281, 145 (1997); 

10. See for example: A. Pich and E. de Rafael, Nucl. Phys. B 58, 311 (1991), S. Bertolini, J.O. Eeg and 
M. Fabbrichesi, Nucl. Phys. B449 (1995) 197, V. Antonelli, S. Bertolini, J.O. Eeg, M. Fabbrichesi 



and E.I. Lashin, Nucl. Phys. B469 (1996) 143, S. Bertolini, J.O. Eeg, M. Fabbrichesi, and E.I. Lashin, 
Nucl.Phys. B514 (1998) 63-92, and references therein. 



